Thin film deposition with time varying temperature 
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We study the effects of time-dependent substrate/film temperature in the deposition of a meso- 
scopically thick film using a statistical model that accounts for diffusion of adsorbed atoms with 
coefficients dependent on an activation energy and temperature. For simplicity, a model where only 
adatoms without lateral neighbors can move is studied by computer simulation, with deposition of 
typically 10 4 atomic layers and temperature changes up to 30i(", near or below the room temper- 
ature range. Linearly varying and exponentially converging temperature cases are considered. If 
the temperature decreases during the growth, the global roughness slowly increases at short times, 
but may show a rapid growth with exponent /3 > 1/2 after ~ 10 3 monolayers. In the exponentially 
converging case, the local roughness shows evidence of anomalous scaling, with a large range of the 
anomalous exponent, and in the linearly decreasing case the local roughness may increase faster 
than a power law. If the temperature increases during the growth, a non-monotonic evolution of the 
global roughness may be observed, with a maximum in the linearly increasing case and a sequence 
of maximum and minimum in the case of exponential convergence to a final temperature. This is 
explained by the competition of kinetic roughening and the smoothing effect of increasing diffusion 
coefficients. If some of these features are observed in experiments, our results suggest to investigate 
the possibility of time varying temperature before relating them to more complex physico-chemical 
mechanisms. 

PACS numbers: 68.55.-a, 68.35.Ct, 81.15.Aa , 05.40.-a 



I. INTRODUCTION 

The morphological properties of a growing film by 
vapor techniques are mainly determined by the bal- 
ance between substrate/film temperature and pressure. 
High temperature helps the system to attain equilibrium 
states, for instance favoring formation of smooth surfaces 
in homoepitaxy, while the increase of pressure leads to 
faster adsorption and drives the system far from the equi- 
librium conditions [l|. For these reasons, many statistical 
models of thin film growth represent the competition of 
the external flux of atoms or molecules and temperature- 
dependent surface processes, such as diffusion, aggrega- 
tion, and reactions d-d]. They account for an external 
flux of atoms or molecules with rate F and for surface 
diffusion of aggregated species with coefficient D, so that 
the ratio R = D/F is the main parameter that repre- 
sents the temperature/pressure interplay (possibly other 
activation energies are also model parameters) [H, 0] ■ 

The main quantity to characterize surface morphology 
is the roughness, which measures thickness fluctuations 
along the film surface. It may be measured for the whole 
surface (global roughness) or inside a box that glides over 
that surface (local roughness). The scaling properties 
with time and size are described by dynamic relations 
Q, which are called normal when local and global fluc- 



tuations scale with the same exponents Q and anomalous 
when local and global exponents arc different 

A question of experimental relevance is the effect of 
changing physico-chemical conditions during film growth. 
For example, a recent experimental work on electrodcpo- 
sition of Prussian Blue films interpreted the growth as a 
process dominated by surface diffusion in which the time- 
increasing adsorption rate is responsible for an anomaly 
in the dynamic scaling relation |ljj. The problem of 
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spontaneous fluid imbibition in a porous medium is an- 
other system that shows anomalous scaling (AS) and is 
modeled as interface growth with time-dependent cou- 
plings fl2j . 

This scenario motivated theoretical works on growth 
models where the microscopic rules of aggregation change 
in time. Ref. [l3[ showed that a sudden change in 
the parameters of the Edwards- Wilkinson (EW) equation 
[h"H of interface growth may be responsible for nontriv- 
ial effects in roughness scaling, such as power-law relax- 
ation to steady states. The possibility of AS was also 
discussed in the context of stochastic growth equations 
with time-dependent couplings and was illustrated for 
the EW equation with time-dependent surface tension 
in Ref. [la ]. It showed good agreement with numeri- 
cal results for models of spontaneous imbibition [l|| [l7| ■ 
In lattice models of film growth, a recent realization of 
AS was shown in competitive models with time- varying 
probabilities and involving uncorrelated deposition jl8j |. 

On the other hand, to our knowledge no work have al- 
ready discussed the effects of time- varying conditions in 
realistic models of thin film deposition. For this reason, 
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this paper is devoted to study a model of deposition and 
diffusion of adsorbed species in which the diffusion coef- 
ficients are affected by a time-dependent substrate/film 
temperature. For simplicity, the model assumes that only 
adatoms in terraces can move, while atoms with lateral 
neighbors are permanently aggregated. Despite its sim- 
plicity, the model shows several nontrivial effects on sur- 
face roughness scaling that may help to understand ex- 
perimental results. In cases of decreasing temperature, 
there is evidence of AS for film thicknesses typical of real 
mesoscopically thick films, with rapidly increasing global 
roughness that suggests nontrivial growth exponents. In 
cases of increasing temperature, a non-monotonic evo- 
lution of the global roughness may be observed, possi- 
bly with a saturation for a wide range of film thickness. 
The main advance from previous works is to show those 
features in realistic conditions of film growth, since the 
model may be viewed an approximate description of low- 
temperature physical vapor deposition. 

The rest of this paper is organized as follows. In Sec. 
II, we will define the model, review basic concepts of dy- 
namic scaling of surface roughness, and discuss scaling 
properties for fixed temperature. In Sec. Ill, we will 
show numerical results for various conditions of decreas- 
ing temperature. In Sec. IV, we will show numerical 
results or various conditions of increasing temperature. 
In Sec. V, we summarize our results and present our 
conclusions. 



II. BASIC DEFINITIONS AND CONCEPTS 

A. Model definition 

The model studied in this work was introduced in Rcf. 
fl9j . It is inspired in the Das Sarma and Tamborenea 
model of molecular beam epitaxy [20j . but changed in 
order to allow a dynamic competition of deposition and 
diffusion, with the possibility of adatoms migrating to 
long distances in terraces. It also resembles the models 
of submonolayer growth with critical island size i = 1 |H , 
although here the deposition of a large number of layers 
is allowed. 

The model is of solid-on-solid type, i. e. there is no 
overhang in the film surface. The deposit has a simple 
cubic lattice structure. The deposition begins with a flat 
substrate of linear size L in the xy plane, with the lat- 
tice parameter being the length unit. The deposit conse- 
quently has Lx L columns and the column height is the z 
coordinate of the topmost adatom. There is an external 
flux of F atoms per site per unit time, measured in num- 
ber of monolayers (ML) per second. An incident atom 
is adsorbed upon landing above a previously deposited 
atom or a substrate site. All adsorbed atoms with no 
lateral and no upper neighbor (i. e. isolated atoms in 
terraces) diffuse with coefficient D. If an adatom has a 
lateral or an upper neighbor, then it is permanently ag- 
gregated at that position. This is a reasonable assump- 



tion for low temperatures and reduces simulation times 
compared with models where all atoms are mobile. It 
has the disadvantage of not respecting detailed balance 
conditions, but various successful applications justify its 
use for submonolayer growth Fig. I illustrates the 
possible steps of some mobile atoms at the film surface. 
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FIG. 1: Illustration of the model: deposition of an atom (1); 
adatom diffusion in terraces (2,4), with possibility of downhill 
movement and irreversible aggregation (3), and with possibil- 
ity of irreversible aggregation at a step edge (5). 

The diffusion coefficient D is expected to depend on 
temperature as D = vq cxp (— E/HbT), where vq is a 
jump frequency and E is the activation energy in a flat 
surface, which typically amounts to tens of eV . The pa- 
rameter R = D/F determines the film properties, since 
this is the average number of steps of a mobile adatom 
during the time of deposition of one layer of atoms (L 2 
atoms). The other relevant parameter is the film thick- 
ness t, proportional to the deposition time. Our model 
allows downhill movement at terrace edges without ad- 
ditional energy barriers (Ehrlich-Schwebel) [2l[. This 
movement always lead to permanent aggregation at the 
step edges. For this reason, the inverse process (uphill 
movement) is not possible. Our simulations are per- 
formed in lattices with L = 512. A fixed deposition 
parameter is vq/F = 10 14 and the chosen activation en- 
ergy value is E = 0.6 cV, which are close to estimates 
for Pd/Pd(100) submonolayer epitaxy 0,[22|. This can- 
not be taken as a direct application due to the different 
lattice structures. On the other hand, the physically im- 
portant parameter is R, which is a function of E/T, thus 
choosing a fixed value of E and changing T does not re- 
strict our conclusions. Since the model is suitable for 
low temperatures, we will consider temperature ranges 
in which R typically varies between 10 1 and 10 5 . Addi- 
tional reasons for this choice are presented in Sec. Ill CI 
In all cases, the temperature is uniform across the film. 

During the simulation of a given sample, the adsorp- 
tion of each new atom takes place in a time interval 
1/ (FL 2 ). After this process, R/L 2 steps of randomly 
chosen free atoms are performed. Since R/L 2 is usually 
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not integer, and it may be small for large L, we keep its 
fractional part for determining the number of steps after 
the next adsorption events. If the temperature changes 
during the film growth, then the value of R is updated 
after the deposition of each complete layer (L 2 atoms), 
which corresponds to a unit time interval. 



B. Dynamic scaling of surface roughness 

The global roughness (or interface width) is defined as 

-2\ V2 



W(L,t) = ({h-h) 



(1) 



where the overbars represent spatial averages and the 
angular brackets represent configurational averages over 
various samples. For short deposition times (small thick- 
nesses), effects of the finite substrate size L are negligible 
and W scales as 



W ~ t , 



(2) 



where /3 is called the growth exponent. For comparison 
with experimental data, the local roughness w [r, t) is 
more useful. It is also defined as a root mean square 
height fluctuation, but the spatial average is limited to a 
box of size r, with r < L, and the configurational average 
is performed by gliding this box over the film surface and 
over various samples. In the regime of thin film growth 
(negligible effects of substrate size) , it follows the Family- 
Vicsek (FV) dynamic scaling relation Q 



w (r, t) « Ar°" g 



Bt 



(3) 



Here, a; is the local roughness exponent (sometimes 
called Hurst exponent), z is the dynamic exponent, B is 
a constant, and g is a scaling function such that g(x) ~ 1 
for x ^ 1 (small box size) and g{x) ~ x~ ai ^ z for x <C 1 
(large box size, where the local roughness coincides with 
the global one). 

In systems with normal scaling, the amplitude A in 
Eq. (J3j) is a constant (time- independent), thus j3 = cti/z 
and the global roughness exponent a equals the local one. 
A scaling analysis of w (r, t) of several growth models in 
one- and two-dimensional substrates is presented in Ref. 
p3j . In systems with AS, the amplitude A scales as 



A~t" 



(4) 



with k characterizing the degree of anomaly [8(. A FV 
relation can also be defined for the global roughness and 
involves the exponent a, but it is useful only if very long 
times / large thicknesses are attained. 



C. Previous works 

Here we restrict the discussion to results on two- 
dimensional substrates (growth in 2 + 1 dimensions), 



which is the subject of the present work. Related full dif- 
fusion models were previously studied numerically with 
fixed temperature during the deposition. These models 
allowed diffusion of atoms with lateral neighbors, with an 
increase in the energy barrier proportional to the num- 
ber of those neighbors. Some authors proposed that they 
had AS [iH or a log arithmic scaling of amplitude A in 
the FV relation [25(. Wilby ct al [26| showed that the 
growth exponent /3 was close to value of the universal- 
ity class of the fourth order nonlinear gro wth equation of 
Villain, Lai and Das Sarma (VLDS) [27|, HI]. Recently, 
renormalization studies explained the long crossover to 
VLDS scaling [H,[IJ. 

The model presented in Sec. Ill Al was recently studied 
numerically, showing a FV relation for t he g lobal rough- 
ness that includes the independence as [lj| 



W = —f(R y t/L z ), 



(5) 



with a « 2/3, z « 10/3, x 0.5, and y « 1. The 
exponents a and z agree with those of the VLDS class [H, 
3lT l32j . The values of exponents x and y were explained 
by scaling arguments j!9| . 

For R > 10 5 , the surface roughness of that model is 
very small even with very large thicknesses (10 or 10 5 
ML) [l9j]. Thus, it is not in the regime of dynamic scal- 
ing. This is an additional reason for our simulations to be 
performed with relatively small values of R (for compar- 
ision, full diffusion models with values of R of the same 
order show larger values of the roughness, particularly if 
ES barriers are included [33|). 



D. Scaling of local roughness 

A recent work on numerical integration of the VLDS 
equation showed evidence of AS in 2 + 1 dimensions [13] , 
which raises the question on a possible AS in our lattice 
model. This feature was not investigated in Ref. [lj|, 
thus it is considered here. 

We performed simulations of the model with 10 < R < 
1000 for fixed R (fixed temperature). Fig. 2 shows the 
local roughness as a function of box size for three values of 
that parameter. For R = 10, there is a split in the curves 
for small r, which suggests AS. However, that split tends 
to disappear as the thickness increases. For instance, 
from t = 4000 M L to t = 8000ML, the change in the local 
roughness for r — 5 is only 7%. For R = 100 and R = 
1000, the split of the curves for different times is almost 
negligible for t > 10 3 . For those reasons, we understand 
that there is no asymptotic AS in our lattice model with 
fixed R, and an apparent anomaly appears only for small 
R and for small thicknesses due to some type of scaling 
correction. 

Fig. 2 also shows that, for small r, the local roughness 
scales with local exponent a; w 1/3 [Eq. {3J|], which 
is significantly below the VLDS value ai = a rj 2/3 
[32j ■ This discrepancy is frequently observed in thin film 
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FIG. 2: Local roughness as a function of box size for R=10 
(full line), R=100 (dashed line) and R=1000 (dashed dotted 
dotted line) at times t = 500ML, 1000ML, 2000ML, 4000ML 
and 8000ML, from bottom to top. The dotted line indicates 
the slope value of 1/3. 



growth models due to scaling corrections, as explained 
in Ref. (23j, and may not be viewed as a failure of the 
predicted model class. 



III. DEPOSITION WITH DECREASING 
TEMPERATURE 



These are situations with slow down of the diffusion 
process, thus roughening is facilitated. First we consider 
linearly decreasing temperature: (I) from Tj = 300A" 
to T F = 270K and (II) from Tj = 27QK to T F = 240K, 
during the time interval for deposition of t max = 10 4 ML, 
which is a reasonable thickness for mcsoscopic films (Tj 
and Tp denote initial and final temperatures, respec- 
tively) . Fig. 3a shows the time evolution of the param- 
eter R and of the temperature, and Fig. 3b shows the 
evolution of the global roughness of the films. For small 
thicknesses (t < 10 3 ML), the slopes of the plots in Fig. 
3b [effective exponents /3; Eq. @] are close to or slightly 
below the VLDS value (3 « 0.2. This may be surpris- 
ing because decreasing temperature favors roughness in- 
crease and was expected to provide a larger (3. However, 
for case (I) , the roughness is very small (typically smaller 
than unity), thus this is not a true scaling region where 
estimates of (3 can be reliably interpreted. For larger 
thicknesses, W rapidly increases in cases (I) and (II). 
For 5 x 10 3 A/L < t < 10 4 ML, the effective exponents 
are (3 ~ 0.70 and (3 0.82, respectively. This nontrivial 
growth exponent cannot be interpreted by existing theo- 
ries of interface growth which assume time-independent 
couplings. 

Recall that (3 = 1/2 characterizes random uncorre- 
cted deposition 0] and is expected when R — > 0. How- 
ever, this not the case here: R ~ 20 at Tp = 24QK 
and R ~ 700 at Tp = 270K. Usual interpretations of 
f3 > 1/2 are preferential aggregation at surface hills in- 
stead of valleys, possibly with some surface instability, 



shadowing effects [35|, |36[, or Ehrlich-Schwebel barriers 
[U . However, these features are also not present in our 
model. Moreover, large (3 is not related to large val- 
ues of the roughness. For instance, very rough deposits 
produced by grain deposition models [371 ] have (3 ~ 0.1, 
which is well below the asymptotic value of the Kardar- 
Parisi- Zhang class [38j . 
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FIG. 3: (a) Parameter R as a function of the thickness, con- 
sidering E = 0.6eV, uq/F = 10 14 , and temperature variation 
of cases (I) and (II). (b) Global roughness as a function of the 
thickness for cases (I) and (II). In (a) and (b), squares cor- 
respond to case (I) and circles to case (II). (c) Local rough- 
ness as a function of box size for cases (I) (solid lines) and 
(II) (dashed lines) at times t = 500ML, 1000ML, 2000ML, 
4000ML and 8000ML, from bottom to top. The inset shows 
log 10 (w) as a function of log 10 (t) for r = 10. 

The exponents j3 > 1/2 obtained here are a special 
consequence of the decreasing deposition temperature. 
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Thus, experimental works which find this result must 
consider the possibility that the temperature is chang- 
ing during the deposition before choosing interpretations 
based on other microscopic mechanisms, such as those 
listed above. The possibility of time-increasing adsorp- 
tion rates must also be considered because it also leads to 
time-decreasing R (this was the case in Ref. [lljj, where 
f3 « 0.5 was obtained). Fig. 3c shows the local rough- 
ness versus box size for some thicknesses in the range 
500ML < t < 8000ML. There is a significant split of 
the curves for small values of r, in contrast to the con- 
stant temperature case shown in Fig. 2 (note that ver- 
tical and horizontal scales are the same in Figs. 2 and 
3c). It suggests AS, with the FV relation [Eq. hav- 
ing time-dependent amplitude A. However, for small and 
constant r, the scaling of A with time does not follow the 
power- law of Eq. §4§ ( scc the inset of Fig. 3c); instead, 
A increases much faster in time. Again, if the slow down 
of surface diffusion is not accounted for, this result may 
naively suggest some sort of surface instability. 
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FIG. 4: (a) Parameter R as a function of the thickness, con- 
sidering E = 0.6eV, vq/F — 10 14 , and temperature variation 
of cases (III) and (IV). (b) Global roughness as a function 
of the thickness for cases (III) and (IV). In all plots, squares 
correspond to case (III) and circles to case (IV). 



Now we consider cases of exponential convergence of 
the temperature to a final value: (III) from Tj — 300A" to 
T F = 270A and (IV) from T 7 = 270A to T F = 240A'. In 
both cases, a characteristic time/thickness t c = 2000ML 



T = Tp + (Tj — Tp) cxp (— i/i c ). 



(6) 



Deposition of t = 10 4 i\/L = 5t c layers is considered, thus 
the final temperatures of the films are very close to T F - 
Cases (III) and (IV) represent systems exchanging heat 
by conduction with a colder reservoir (e. g. an initially 
heated system separated from the surroundings by a con- 
ducting material). Fig. 4a shows the time evolution of 
the parameter R and temperature in cases (III) and (IV). 
It contrasts with the linear decay in Fig. 3a because 
R rapidly decreases up to t ~ t c , but decreases slowly 
during most of the deposition time, when the system at- 
tained a temperature close to the surroundings. Fig. 4b 
shows the evolution of the global roughness for cases (III) 
and (IV) . Again, W increases slowly for small thicknesses 
and large effective exponents j3 are obtained for large 
thicknesses. In case (IV), where R decreases almost two 
orders of magnitude during the film growth, an exponent 
(3 > 1/2 is also found, which reinforces the conclusion 
that a careful interpretation of this result in real systems 
is necessary. 
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FIG. 5: Surface topography of films with thicknesses t = 
500ML and t = 8000ML in (top) case (IV) and (bottom) 
with fixed temperature T — 270K. 

Fig. 5 shows the surface topography of films with 
thicknesses t = 500ML and t = 8000MX in case (IV) 
and with fixed temperature T = 270A". While the films 
grown in fixed temperature are always smooth, a remark- 
able roughening is observed in the films grown with de- 
creasing temperature, consistently with the large values 
of /3. 

Figs. 6a and 6b show the local roughness as a function 
of box size for cases (III) and (IV), respectively. Both 
plots show evidence of AS due to the split of the curves 
for small r. Surprisingly, this trend persists for thick- 
nesses t > 2t c , where the temperature is changing very 
slowly. Fig. 6c shows the scaling of w(r = 10, t), which 
gives the anomaly exponents k = 0.28 ± 0.01 for case 
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models in Ref. although those models are not so 

close to real film growth. Finally, we observe that the 
local roughness exponents measured in Figs. 3c, 5a, and 
5b are close to 1/3, similarly to the constant temperature 
case (Sec. Ill PI) . 



IV. DEPOSITION WITH INCREASING 
TEMPERATURE 

First we consider two cases of linearly increasing tem- 
perature: (V) from T[ = 240A" to T F = 270A and (VI) 
from Tj = 270A' to T F = 300A", both from t = to 
tmax = 10 ML. Fig. 7a shows the time evolution of the 
parameter R and of the temperature and Fig. 7b shows 
the evolution of the global roughness. The nontrivial re- 
sult is the nonmonotonic variation of the roughness. It 
increases for short times due to spread of correlations and 
excitation of modes of increasing wavelength, as usual in 
kinetic roughening [H, |39[. However, it attains a maxi- 
mum and begins to decrease due to the smoothing effect 
of large diffusion coefficients. The effect is more pro- 
nounced for smaller temperatures because the relative 
increase of R is larger. 
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FIG. 6: (a), (b) Local roughness as a function of box size 
for cases (III) and (IV), respectively, at times t = 500ML, 
1000ML, 2000ML, 4000ML and 8000ML, from bottom to top. 
(c) log 10 (uj) as a function of log 10 (t) for r — 10 in cases (III) 
(squares) and (IV) (circles). The dashed lines indicate the 
linear fits for the k evaluation. The corresponding values of 
k are also shown. 



(Ill) and k = 0.47 ± 0.02 for case (IV). This is AS found 
in conditions plausible in experimental work (note that 
maximal temperature changes arc 30 A). It is important 
to stress that the AS found in cases (III) and (IV) is an 
apparent scaling restricted to the thickness range studied 
here. At long times, the deposit will attain the character- 
istic features of films grown with fixed T F , which do not 
have AS, but have normal VLDS scaling. Cases (I) and 
(II), with linearly decreasing temperature, are necessarily 
limited to a finite time range, thus there is no asymptotic 
scaling. These features contrast with the true AS of the 
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FIG. 7: (a) Parameter R as a function of the thickness, con- 
sidering E = 0.6eV, vo/F = 10 14 , and temperature variation 
of cases (V) and (VI). (b) Global roughness as a function 
of the thickness for cases (V) and (VI). In all plots, circles 
correspond to case (V) and squares to case (VI). 

The situation with exponential convergence of tern- 
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perature [Eq. (|6]) with Tj < Tp] is more interesting, 
although the variations in W are smaller. We con- 
sider cases (VII) from T/ = 240X to T F = 270K and 
(VIII) from T/ = 270A" to Tp = 300K, again from 
t = to t mal = 10 4 AfL, with characteristic decay time 
i c = 2000ML. Fig. 8a shows the time evolution of the 
parameter R and of the temperature. Fig. 8b shows 
the evolution of the global roughness, which presents 
a maximum at small thicknesses for case (VII) and a 
plateau for case (VIII) . These effects are similar to those 
of the linearly increasing T [(V) and (VI)]. However, for 
larger thicknesses, the roughness slowly increases. This 
is expected because the temperature slowly varies for 
t > t c , thus the usual features of kinetic roughening (c. 
g. roughness increasing with time) are observed. 
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FIG. 8: (a) Parameter R as a function of the thickness, con- 
sidering E = 0.6eV, uo/F = 10 14 , and temperature variation 
of cases (VII) and (VIII). (b) Global roughness as a function 
of the thickness for cases (VII) and (VIII). The arrows in- 
dicate a maximum and a minimum of global roughness. In 
all plots, circles correspond to case (VII) and squares to case 
(VIII). 

The presence of two extremal values in the roughness 
(a maximum and a minimum) was already shown in an 
electrodeposition model where the time-dependent ad- 
sorption rate had a maximum [40j . This corresponds to a 
minimum of an effective diffusion-to-deposition rate R for 
that model. The increase of R after that minimum seems 
to be responsible for the effect of reducing the roughness 
in a certain time window, similarly to our model. The 



deposition with large R may also lead to roughness oscil- 
lations at short times due to fluctuations between partial 
and complete filling of the first atomic layers. This fea- 
ture was shown by Pal and Landau [HI, S3 m some lattice 
models and has no relation with the roughness evolution 
presented here. Indeed, those oscillations disappear af- 
ter deposition of approximately 10 layers (maybe less), 
in contrast to our models. 



V. CONCLUSIONS 

We studied a thin film growth model in which the sur- 
face diffusion coefficients of adsorbed species are related 
to the substrate/film temperature in conditions where 
that temperature increases or decreases during the depo- 
sition. The model assumes that only adatoms in terraces 
can move and was formerly studied only in the case of 
constant deposition temperature. 

In all cases analyzed here, 10 4 atomic layers were grown 
with temperatures in the range of room temperature 
or below and variations during the growth up to 30K, 
which arc feasible conditions with many experimental 
techniques. Cases of temperature varying linearly in time 
and of an exponential convergence to a final temperature 
were separately analyzed, the later representing systems 
exchanging heat by conduction with a thermal reser- 
voir. If the temperature decreases during the growth, 
the global roughness slowly increases at short times but 
eventually turns to be a rapidly increasing function of 
time. In several cases, growth exponents /3 > 1/2 are 
obtained in the thickness range 10 3 ML < t < 10 4 ML. 
This shows that the simple physical mechanism of reduc- 
ing the temperature during the growth may lead to non- 
trivial exponents which are usually related to more com- 
plicated features, such as step energy barriers or other 
mechanisms that generate instabilities. Moreover, the 
local roughness shows evidence of anomalous scaling in 
the thickness range analyzed. In the cases of exponen- 
tial convergence to a final temperature, several values 
of the anomaly exponent may be obtained in the thick- 
ness range analyzed. In the cases of linear temperature 
decrease, the local roughness in a small box size may 
increase in time even faster than a power law. 

If the temperature increases during the growth, a non- 
monotonic evolution of the global roughness may be ob- 
served. In cases of linear temperature increase, a maxi- 
mum of roughness is usually observed, being more pro- 
nounced for lower temperature ranges. In the case of 
exponentially converging temperature, a maximum and 
a minimum may be observed (eventually turning into a 
plateau in the logW x logi plot). These features are 
interpreted as a consequence of competition of kinetic 
roughening, which leads to roughness increase, and the 
smoothing effect of increasing diffusion coefficients. 

Stochastic growth equation approaches [ll| and lattice 
models [l8| have already considered the effect of time- 
varying couplings in interface growth, with emphasis on 
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anomalous scaling properties. The main advance of the 
present work is to show that a series of nontrivial fea- 
tures of surface roughness scaling (including anomalous 
scaling) appear in realistic conditions of film growth. We 
expect that our results may be helpful for the interpreta- 
tion of previous or future works on thin film deposition 
in which the temperature may vary during that process. 
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